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1.0 SUMMARY 


The program L217 (EOM) is structured as a system of overlays, one main overlay with 
four primary overlays and their associated secondary overlays. 

Input into the program is made via cards and magnetic files (tapes or disks). Output 
from the program consists of printed results and magnetic files containing matrices of 
equations of motion coeflicients and aerodynamic forces. 

Although L217 (EOM) serves as a module of the DYLOFLEX system, it can be operated 
as a standalone program. Subroutines used by L217 (EOM) include routines obtained 
from the standard FORTRAN library, and routines from the DYLOFLEX alternate 
subroutine library. 


2,0 INTRODUCTION 


The computer program L217 (EOM) was developed for use as either : a standalone 
program or as a module of a program system called DYLOFLEX (ref. 1) (see fig. 1). The 
Equations of Motion program (L217) was designed to meet the DYLOFtLEX contract 
requirements as defined in reference 2. These requirements specify the need for a 
program that can assemble and generate data needed to formulate "the equations of 
motion for a flight vehicle, and the necessary aerodynamic forces that can be used in 
the flight loads analysis of that vehicle. \ \ \ 

The program was designed, coded, and documented according to the ^programming 
specifications developed for DYLOFLEX^. , ; . . 

The objective of this volume is to aid those persons who will maintain and/or modify the 
program in the future. To meet this objective the following items are discussed: 

• Program design and structure ^ 

• Purpose and operation of each overlay . cf J 

• Input, output, and internal data bases / . 

• Test cases 


Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads System to 


Interface With FLEXSTAB. Boeing-NASA Contract NASl-13918. BCS-G0701, September 1975, 
(internal report) 
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Figure 1. - DYLOFLEX Flow Chart 
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3.0 PROGRAM DESIGN AND STRUCTURE 


The program is structured as a system of overlays, the main overlay and four primary 
overlays with their associated secondary overlays (see’ fig. 2). 


Main overlay : 

(L217,0,0) 

Program L217vc 

Primary overlay 

(L217,2,0) 

Program STRUCT 

Primary overlay 

(L217.3.0) 

Program FLXAIC 

Primary overlay 

(L217,4,0) 

Program DLAIC 

Primary overlay 

(L217,5,0) 

Program DLPRES 


The main overlay named L217vc (the v and c are version and correction identifiers. 
See sec. 3.1) reads' data cards that direct the execution of the primary overlays. It also 
aids communication between the primary overlays via labeled common blocks. 

The 2,0 primary overlay, STRUCT, generates the equations of motion matrices that 
are dependent upon the structural data only. They are: 

[Mx] Generalized stiffness matrix 
[M2] Generalized damping matrix 
[M3] Generalized inertia matrix 

The remaining primary overlays generate the equations of motion matrices that are 
dependent upon the aerodynamic data. They are: 

[M4] Generalized aerodynamic stiffness matrix 
[M5] Generalized aerodynamic damping matrix 

{(l) The distance(s) from the gust reference point to the control point(s) of the 
gradual penetration gust zones. 

[^] The generalized gust forcing function matrix 

The 3,0 primary overlay, FLXAIC, generates the aerodynamic equations of motion 
that are based upon the FLEXSTAB (ref. 3) AIC (aerodynamic influence coefficient) 
matrices. 

The 4,0 primary overlay, DLAIC, generates the aerodynamic equations of motion 
based upon Doublet Lattice (L216) (ref. 4) AIC matrices. 

The 5,0 primary overlay, DLPRES, generates the aerod3mamic equations of motion 
based upon generalized forces and unsteady pressures from Doublet Lattice (ref. 4). 




L21 7.5,6 
DLPEOM 


L21 7,5,5 
DLPODR 


1217.4,11 

DLAEOM 




























Although L217 (EOM) serves as a module of the DYLOFLEX system, it is operated as 
a standalone program. When the program is run by itself, it becomes the user’s 
responsibility to generate input data in the format required by the program. (See vol. 
I of this document for a description of the required data and formats.) L217 (EOM) is 
heavily dependent upon data generated prior to its execution. Figure 3 displays the 
different files, which L217 (EOM) may require, according to user specified options. The 
following programs produce or read the files displayed in figure 3. 


Program Reference number 

3 

4 

5 

6 

7 

8 

The program requires subroutines which are not part of the L217 (EOM) code. Some 
routines are automatically obtained from the standard GDC FORTRAN library when 
the program is loaded. Others, however, are stored in the DYLOFLEX alternate 
subroutine library that must be declared at the time of loading. Subsequent sections 
describe each overlay separately and contain tables displaying the routines called and 
the library in which they are located. All subroutines have one entry point. 

This volume describes the program in a macro sense. A more detailed discussion 
appears in the comments contained in the program source code. Each routine contains 
a preface describing the routine’s purpose, author, analytical steps, modification 
history, input data, output data, glossary of variables, and a list of other routines 
called. Embedded within the executable code are comments labelling each section and 
explaining logical branches. 


FLEXSTAB 

L216(DUBFLX) 

L215(INTERP) 

L22KTEV156) 

L219(EQMOD) 

L218(LOADS) 


3.1 OVERLAY(L217,0,0), L217VC 
The main overlay of L217 (EOM) is named L217vc where: 

V Is a letter indicating the program version 

c Is an integer indicating the correction set applied to the v version 

Purpose of L217vc 

L217vc performs certain bookkeeping tasks, directs the execution of the primary 
overlays, and aids communication between primary overlays via labeled common 
blocks. 



FLEXSTAB 


SDINDX 


SDDATA 


From 

L216(DUBFLX) 


NAETP 


From 

L215 0NTERP) 


GOT APE 

NGETP 

GSTIFF 

GMASS 

SATAP 

L .. . ■ 

... . 

L V 

J 

V 


L217 (EOM) 


Printed 

output 


EOMTAP , 

To L221 (TEV156) 
or L219(EQMOD) 


EOMLOD 
To L218(LOADS) 


Figure 3. -- L21 7 (EOM) Communication Via External Files 





Analytical Steps of L 217 vc 


The L217vc performs its tasks in the following steps: 

1. The subroutine PRGBEG is called to place the program header on the printed 
output. 

2. A data card is read. It must begin with $EQM to ensure that the input file is 
correctly positioned. If it does not contain $EOM, execution is terminated. 

3. A program directive card is read, printed, interpreted, and acted upon according 
to the following logic: 

a. If the keyword is $TITLE jump to step 3 again. 

b. If the keyword is $NONSTOP set the nonstop switch and jump to step 3 
again. 

c. If the keyword is ^STRUCTURAL, jump to step 4. 

d. If the keyword is $AERO followed by FLEXSTAB jump to step 5. 

e. If the keyword is $AERO followed by DOUB...A1C jump to step 6. 

f. If the keyword is $AERO followed by DOUB...PRES jump to step 7. 

g. If the keyword is $QUIT jump to step 8. 

4. Overlay (L217,2,0), STRUCT is called. When it is finished; control returns to 
step 3. 

5. Overlay (L217,3,0), FLXAIC is called. When it is finished, control returns to step 

3. . 

6. bverlay (L217,4,0), DLAIC, is called. When it is finished, control returns to step 
3. 

7. Overlay (L217,5,0), DLPRES, is called. When it is finished, control returns to 
step 3. 

8. Subroutine PRGEND is called to place the program termination message on the 
printed output file, and the execution is terminated. 

If errors were detected, the program jumps to step 8. 

The macro flow chart of this overlay is shown in figure 4. 

The routines called are displayed in tidile 1. 




Figure 4. — Macro Flow Chart of Overlay (L217, 0, 0), L217 vc 



















Table 1. — Routines Called By L217vc 

OVERLAY (L2 17, 0,0) 

PFOGFAM L217VC 


DLAIC 

OVERLAY (L2 17, 4,0) 

DLRRES : - . 

; OVERLAY (L2 17, 5,0) 

FLicAIC 

OVERLAY (L2 17, 3,0) 

FNDKEY+ 


IPDCRD+ 

■■ - 

PFGBEG* 


PRGEND* 


RETURNF+ 


STRUCT 

OVERLAY<L217,2,0) 


+ indica-tes a routine obtained from 

the DYLOFLEX alternate subroutine library* 


All Others are local to L217(EON). 
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I/O Devices of L217vc 


The L217vc reads program directive cards and. writes them on the printed ouput file. 
All other I/O accomplished by L217 (EOM) is done in lower level overlays. 

3.2 OVERLAY(L217,2,0), STRUCT 

Purpose of STRUCT 

The L217 (EOM)’s first primary overlay is named STRUCT. STRUCT generates the 
equations of motion matrices [Mj], [^ 2 ], and [M 3 ]. 

Analytical Steps of STRUCT 

STRUCT performs its task in the following steps: 

1 . Initialize the error count and print a message indicating that STRUCT is being 
executed. 

2. Call overlay (L217,2,l) RSTRCT to read the card input data directing the 
structural data processing. 

3. If control surface freedoms are to be added to the existing generalized mass and 
stiffness matrices, call overlay (L217,2,2), FORMAC, to form the matrices [A] 
and [C]. (See sec. 3.2.2.) 

4. Call overlay (L217,2,3), FORM13, to generate the matrices [Mj], [l 8 l 2 ]> 

[M 3 ] and save them on the temporary storage file named S2L217. 

5. Print a message indicating that STRUCT is finished and return to the main 
overlay, L217. 

If errors were detected in the secondary overlays, the program jumps to step 5. 

The macro flow chart of this overlay is shown in figure 5. 

The routines called are displayed in table 2 . 

I/ODevices of STRUCT 

STRUCT prints two messages on the output file; NOW IN STRUCT.. .and STRUCT IS 
FINISHED... .All other file communication in STRUCT is performed by its secondary 
overlays. Figure 6 shows the files providing input and output to STRUCT’s overlays. 

In subsections 3.2.1 through 3.2.3, the matrices read from and written onto each file 
are listed. (See sec. 6.4 of vol. I, Engineering and Usage ^ for detailed description of the 
format and content of files used for external communication.) Section 3.6.3 of this 
volume contains the same information for files used only within L217 (EOM) for 
temporary storage. 
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ENTRY 



{ Read structural 
card input data. 


( Form [A] and [C] 
only if control 
surface freedoms 
are being added. 


I 


Form [M.^l - [Mg] 


Figure 5. — Macro Flow Chart of Overlay (L217, 2, 0), STRUCT 




Table 2. — Routines Called by STRUCT 

6VEFLAY(L217,2,0) 

PFCX3RAM STFUCT 


FORMAC - OVERLAY(L217,2,2) 
FORM 13 - OVERLAY (L2 17, 2, 3) 
RSTRCT - OVERLAY (L2 17, 2,1) 



L217, 2,0 
STRUCT 



Figure 6, — Files Providing 


] 



[M*2l (M3] 


to STRUCT Overlays 



3.2.1 OVERLAY(L217,2,l), RSTRCT 

Pufpose of RSTRCT 

RSTRCT reads the card input data directing the processing of structural data. Options 
and constants are stored in labelled common blocks for later use. If either the 
generalized mass or generalized stiffness matrix is read from cards, it is written onto 
the temporary storage file named S1L217 for later retrieval by overlay (L217,2,3), 
FORM13. 


The macro flow chart of this overlay is shown in figure 7. 

The routines called are displayed in table 3. 

I/O Devices of RSTRCT 

RSTRCT reads card images from the input file and prints them on the output file 
along with descriptive comments and any necessary diagnostics. 

When either the generalized mass matrix or generalized stiffness matrix is input on 
cards, RSTRCT writes them onto the temporary storage file named S1L217 (see fig. 6). 

3.2.2 OVERLAY(L217,2,2), FORMAC 

Purpose of FORMAC 

FORMAC is called to generate the matrices [A] and[C]. The two matrices will 
become partitions of [Mg], the generalized mass matrix. They represent 
control-surface, generalized coordinates being added to the original structural data. 

Anal 3 rtical Steps of FORMAC 

FORMAC generates [A] and [C] in the following steps: 

1. Open the files ICSTPE (containing the control surface data) and SATAP (the file 
of interpolation arrays). Set aside storage for the arrays [^], [m£s], and 
[CSDATA]. 

Repeat steps 2 through 5 for each control surface adding freedoms; ICS = 1, 
NCSEF. 

2. Read the nodal data for control- surface IS from S1L217. 

3. Call BLDMCS to add the control surface mass and inertia data to [m6s] . 

4. Read the interpolation array {SA} for the control surface from SATAP. 

5. Interpolate for the modal deflection and/or slopes as requested and build the 
matrix [^]. 



ENTRY 


1 
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Table 3. — Routines Called by RSTRCT 


OVERLAY (L217, 2, 1) 
PFOGFAM RSTRCT 


DELETR+ 

FETAD+ 

FETDEL+ 

INITIR+ 

IFDCRD+ 

MATCFD 

NAMFIL+ 

F DCS CD 

SETMFL 

STAFTR+ 

WFTETP+ 


IRDCRD+ 


+ indicates a routine from the DYLOFLEX 


alternate subroutine library. 



6 . Close the flies ICSTPE and SATAP. 

7. Initialize the file S1L217 and generate the temporary array. 

[PCSTM] = [0csl^ 

8 . Generate [A] and write it on S1L217. 

[A] = [PCSTM] [0st1* 

9. Generate [C] and write it on S1L217. 

[C] = [PCSTM] [0csl. ; 

10. Return to the calling program. 

The macro flow chart of this overlay is shown in flgure 8 . 

The routines called are displayed in table 4. 

I/O Devices of FORMAC 

FORMAC’s communication via external files is displayed in flgure 6 . 

FORMAC reads from the flies: 

S1L217 the control surface nodal data written there by RSTRCT 

SATAP the modal interpolaton array(s) specifled by surface number on card input 

FORMAC writes onto the file: 

S1L217 the matrices [A] and [C] 

3.2.3 0VERLAY(L217,2,3), FORM13 

Purpose of FORM 13 

FORM13 generates the structural equations of motion ([Mi], [M 2 ], and [M 3 ]) and 
saves them on the temporary storage file named S2L217. 

Analytical Steps of FORM13 

FORM 13 performs its tasks in the following steps: 

1 . Initialize the output file S2L217 and variably dimensioned array storage. 
Initialize the array [M 2 ], if structural damping was specifled by the user. 
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ENTRY 


Initialize the files 
ICSTPE and SATAP 
and the arrays 
[0] [l/CS] [CSDATAl 


ICS-1.NCSEF 


Loop on control surfaces 
adding new generalized 
coordinates. 


S1L217 


Call READTP to 
read [CSDATA] . 


BLDMCS to o 
construct [MCS] 


ICall READTP to 
SATAP ^SA). 


'Coordinates^ 


Reference 



Figure 8. — Macro Flow Chart of Overlay (L217,2,2), FORMAC 


20 





Table 4. — Routines Called by FORM AC 


OVERLAY (L217,2,2) 
PROGRAM PORMAC 


AINTG+ 

AINTL + 

BLDMCS 

GATED 

DELETR+ 

FETAD+ 

FETDEL+ 

INITIR+ 

MATMLT 

PRINTM+ 

READTP+ 

SETMFL 

STARTR+ 

WRTETP+ 


+ Indicates a routine from the DYLOFLEX 


alternate suroutine library 



2. Initialize [Mj], read the generalized stifTness matrix into its upper left partition, 
scale if requested, and write[Mi] onto the file S2L217. 

3. If structural damping was requested, read the damping factors into [ 1 ^ 12 ] and 
multiply them by the square root of the corresponding diagonal elements of 
[Ml]. 

4. Erase [Mj], initialize [M 3 ], and read the generalized mass matrix into the 
upper left partition. 

5. If structural damping was requested, multiply the elements of [lS[ 2 l by the 
square root of the corresponding diagonal elements of [M 3 ], and wrile [M 2 ] onto 
S2L217. 


6 . Read the matrices [A] and [C] from S1L217 and store in [M 3 ] as shown in the 
following equation, and write [M 3 ] onto S2L217. 


[M 3 I 


[M][A]T 
[A] (Cl 


Note: [a] and [C] are present only when new generalized coordinates are being 
added. 


7. Write an end of file on S2L217, close the file, and return to the calling program. 
The macro flow chart of this overlay is shown in figure 9. 

The routines called are displayed in table 5. 


I/O Devices of FORM 13 

FORMlS’s communication via files is displayed in figure 6 . 

FORM13 reads from the flies: 

GMASS Generalized mass matrix (if not read from cards). 

GSTIFF Generalized stiffness matrix (if not read from cards). 

S1L217 Matrices [A] and [C] (if control surface freedoms are being added) and the 
matrices[M] and [K] if they were read from cards by RSTRCT. 

FORM13 writes onto the file: 

S2L217 the matrices [Mi], [M 2 ], imd[M 3 ] 
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Figure 9. — Macro Flow Chart of Overlay (L217, 2, 3), F0RM13 
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S2L217 


FETDEL 
closes S2L217 

TXT 

C end 3 


fcall WRTETP to 
Iwrite [M 3 ] . 


Figure 9. — (Concluded) 




Table 5. — Routines Called by FORM 13 


OVEFLAY (L217,2,3) 
PRCX3RAM FORM13 


DELETR+ 

FETAD+ 

FETDEL+ 

INITIR+ 

MATAPE 


MLDIAG 

PRINTM+ 

SCLMLT+ 

SETMPL 

STARTR + 

TRANAE 

WFTETP+ 


^ FETAD+ 
FETDEL+ 
FEADTP+ 


indicates a routine from the DYLOFLEX 


alternate subroutine library 


3.3 OVERLAY(L217,3,0), FLXAIC 


Purpose of FLXAIC 

FLXAIC generates aerodynamic force matrices and equations of motion coefUcient 
matrices based upon the AIC matrices (Ap^) generated by FLEXSTAB (ref. 3). 

Analytical Steps for FLXAIC 

FLXAIC performs its task in the following steps: 

1. Initialize the error count and print a message indicating that FLXAIC is being 
executed. 

2. Call overlay (L217,3,l), FLXCDS, to read FLXAIC card input data. 

3. Call overlay (L217,3,2), FLXGEO, to read geometry data from the FLEXSTAB file 
known as GDTAPE. 


4. Call overlay (L217,3,3), FLXMOD, to interpolate for the modal data required by 
FLXAIC. 

5. Call overlay (L217,3,4), FLXFRC, to calculate the aerodynamic forces at all force 
points. 

6. Call overlay (L217,3,5), FLXODR, to sort the aerodynamic forces and prepare the 
output file known as EOMLOD. 

7. Call overlay (L217,3,6), FLXEOM, to generate the aerodynamic equations of 
motion matrices and write them onto the file EOMTAP after the structural 
equations of motion matrices. 

8. Print a message indicating that FLXAIC is finished and return to the main 
overlay, L217. 

If errors were detected in the secondary overlays, the program jumps to step 8. 

The macro fiow chart of this overlay is shown in figure 10. 

The routines called are displayed in table 6. 

The equations solved by FLXAIC are summarized as follows: 
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Figure 10. —Macro Flow Chart of Overlay (L21 7,3,0), FLXAtC 



Table 6. — Routines Called by FLXAIC 

OVEFLAY(L217,3,0) 

PFOGRAM FLXAIC 


FLXCDS 

— 

OVEFLAY (L217,3,1) 

FLXGEO 

- 

OVERLAY (L2 17, 3, 2) 

FLXMOD 

- 

OVERLAY (L217,3,3) 

FLXFPC 

- 

OVERLAY (L217,3,4) 

FLXODR 

- 

OVERLAY (L217,3,5) 

FLXEOM 

- 

OVERLAY (L217,3,6) 


FETAD+ 

FETDEL+ 


+ indicates a routine from the DYLOFLEX 


alternate subroutine library 



[Fj] =-q [FAC][2][A][AIC][«'cp] 

[F 2 ] =4[FAC][2][A][AIC][0CPl 

[FG] = q [FAC] [2] [A] [AIC] [RG] [GC] 

[M4] =-[^Fpl^tSF][Fil 
[M5] =-[0Fpl’'^tSF][F2l 
[?] = [0FP1’^ISF][FG] 

An * indicates that the matrix product is not performed in the normal fashion. 


For instance, [2] is actually stored as NTOTF x four but in the matrix multiply, it acts 
like a matrix of size NAFP x NTOTF. The [GC] is stored as NACP x one but it acts like 
a matrix NACP x NGCP in the matrix multiply. (See description of the file S1L217 just 
ahead of Hg. 36.) 

I/O Devices of FLXAIC 

FLXAIC prints two messages on the output file; NOW IN FLXAIC. ..and FLXAIC IS 
FINISHED... .All other file communication in FLXAIC is performed by its secondary 
overlays. Figure 11 shows the files used in FLXAIC’s secondary overlays. 

In subsections 3.3.1 through 3.3.6 the matrices read from and written onto each file are 
listed. (See secs. 6.4 to 6.5 of vol. I, Engineering and Usage, for detailed descriptions of 
the format and content of the files used for external communication.) Section 3.6.3 of 
this volume contains the sa,me information for L217 (EOM)’s temporary scratch files. 


3.3.1 OVERLAY(L217,3,l), FLXCDS 

Purpose of FLXCDS 

FLXCDS reads the card input data directing the L217 (EOM) generation of equations of 
motion matrices from FLEXSTAB AIC matrices. (See card set 3.0 in vol. I of this 
document. Engineering and Usage.) Options, constants, and file names are stored in 
labelled common blocks for later use. The card data defining the correspondence 
between FLEXSTAB bodies (slender bodies and thin bodies), and the L215 (INTERP) 
modal interpolation surfaces are stored in the arrays {BODTAP} and {ICSTAB} and 
written onto the temporary storage file S1L217. 

The macro flow chart of this overlay is shown in figure 12. 

The routines called are displayed in table 7. 
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Figure 11. — Files 


1217 , 3.0 

FLXAIC 



Input/Output to FLXAIC Overlays 






















Table 7. — Routines Called by FLXCDS 


OVEBLAY (L217,3,1) 
PROGRAM FLXCDS 


FNDKEY* 

IRDCRD+ 

IFOL+ 

NAMFIL+ 

WFTETP+ 


+ indicates a routine in the DYLOFLEX 


alternate subroutine library 



I/O Devices Used 


FLXCDS reads card set 3 from the input file and prints the card images on the output 
file with descriptive comments and any necessary diagnostics. 

FLXCDS writes onto the file: 

S1L217 Arrays {IBDTAB} and {ICSTAB} (see sec. 3.6.3 for contents) 

3.2.2 OVERLAY(L217,3,2), FLXGEO 

Purpose of FLXGEO 

FLXGEO reads the FLEXSTAB geometry data from the file CIDTAPE, combines it with 
the card input data saved on S1L217, and prepares the geometry data required 
thoughout FLXAIC. The matrices produced are written onto the file S1L217 
(overwriting the arrays of card input data). See section 3.6.2 for a description of the 
file’s contents and format. 

The macro flow chart of this overlay is shown in figure 13. 

The routines called are displayed in table 8. 

I/O Devices of FLXGEO 


FLXGEO reads from the files: 

S1L217 Arrays {BODTAB} and {ICSTAB} written there by FLXCDS 
GDTAPE FLEXSTAB geometry data (see sec. 6.4.6 of vol. I) 

FLXGEO writes onto the file: 

S1L217 Geometry matrices and arrays directing the FLXAIC processing (see 
sec. 3.6.3 for content and format) 

3.3.3 OVERLAY(L217,3,3), FLXMOD 

Purpose of FLXMOD 

FLXMOD uses the SA arrays read from the file SATAP (produced by L215 (INTERP)) to 
find the following matrices: 

[i^Cp] Modal deflections at the aero control points (thin body box quarter chord 
and slender body segment midpoints) 


[‘^►'CP] 

[«^FP] 


Modal slopes at the aero control points 
Modal deflections at the aero force points 


ENTRY 


Initialize geometry 


RFLXGE reads the 
FLEXSTAB geometry 
from GDTAPE 












Figure 13. — (Concluded) 



Table 8, — Routines Called by FLXGEO 

OVERLAY (L217, 3,2) 

PFOGRAK FLXGEO 


GTOLT* 

PRTGCP 

PRTGFP 

RFLXCD 

RFLXGE 


SUMSIB 

WFTGEO 


< IPOL+ 
READTP+ 
^ FETAD+ 
FETDEL* 
FSF+ 
GZONES 

WRTETP+ 


* indicates a routine in the DYLOFLEX 


alternate subroutine library 



The DYLOFLEX library routine AINTL is called to interpolate for the modal data using 
the SA array(s) from SATAP and the nodal coodinates read from S1L217. 

The macro flow chart of this overlay is in figure 14. 

The routines called are displayed in table 9. 

I/O Devices of FLXMOD 


FLXMOD reads from the files: 

S1L217 The matrices 
{SURFTB} 

{IPNSRF} 

{X} 

{Y} 

{IPNSCP} 

{XCP} 

{YCP} 

SATAP SA arrays from L215 (INTERP). (see vol. I of this document and ref. 5 for a 
description) 

FLXMOD writes onto the file: 

S1L217 Matrices [<^cp]» [‘^’cp]» [<^*FP] 

Note that if the number of modes is greater than 20 the matrices are partitioned and 
the trio is repeated as often as necessary to include all modes. 

8.3.4 OVERLAY(L217,3,4), PLXFRC 

Purpose of FLXFRC 

FLXFRC generates the aerodynamic forces matrices [Fi], [F 2 ], and [FG] where: 

[Fj I = -q [FAC] [2] [A] [AIC] [0'cp] 

[F 2 ] =4[FAC)[2][A][AlC][0cpl 

[FG] = q [FAC] [2] [A] [AIC] [RG] [GC] 

If the number of modes exceeds 20, the matrices [^cp] [<^'cp] read in 

partitions, and the matrices [Fi] and [F 2 ] are generated in partitions with the same 
number of columns. 

The * above the matrices [2] and [GC] indicates that the operation performed is not a 
normal matrix multiply. Both matrices are described in section 3.6.3 with the file 
S1L217. 


all are described in 
section 3.6.3 
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. Figure 14. — Macro Flow Chart of Overlay (L21 7, 3, 3), FLXMOD 





Figure 14. — (Concluded) 















Table 9. — Routines Called by FLXMOD 


OVEFLAY (L217,3,3) 
PEOGFAM FLXMOD 


AINTL+ 

BINSRT 

FETAD+ 

FETDEL* 

INITIR+ 

PFINTM+ 

READTP+ 

SETMFL 

STAFTF+ 

WETETP+ 


♦ indicates 
alternate 


a routine in the DYLOFLEX 
subroutine library. 



: V 


The macro flow chart of this overlay is shown in figure 15. 

The routines called are displayed in table 10. 

I/O Devices of FLXFRC 

FLXFRC reads from the files: 

S1L217 Matrices [A], [2], [FAC], [GC], [R*G], [<^cp]- [^'cp]> and [d>Fp] 

(See sec. 3.6.3 for a description.) 

SDINDX Table describing the FLEXSTAB aerodynamic data on SDDATA 
SDDATA Matrix AIC by rows 

The AIC matrix read is known inside FLEXSTAB as Ap^ . (See sec. 6.4 of vol. I and 
ref. 3 for a complete description of both SDINDX and SDDATA.) 

FLXFRC writes onto the file: 

S3L217 the matrices 

[Fi] (partitioned if necessary) 

[F 2 ] (partitioned if necessary) 

[FG] (by columns) 

3.3.5 OVERLAY(L217,3,5), FLXODR 

Purpose of FLXODR 

FLXODR sorts the aerodynamic force matrices according to the interpolation surface 
number corresponding with each element. The resulting surface force matrices are 
written onto the file EOMLOD. 

The macro flow chart of this overlay is shown in figure 16. 

The routines called are displayed in table 11. 

I/O Devices of FLXODR 


FLXODR reads from the files: 


S1L217 The matrices 
{SURFTB} 
{IPNSRF} 
{X} 

{Y} 

{AREA} 


See section 3.6.3 
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Rewind S1 L21 7 and call READTP 

to read {A}, [S],{FAC}, {GC), and {RG} 

from the 1st logical file of SI L217. 


( Loop on partitions 
of modal matrices. 

Call READTP to read [0Qp] |p and [0^Qp] 
from the 2nd logical file of S1L217, 

Rewind SDDATA and 
position SDDATA before AlC. 


I Read row I of [AlC] . 


{ Call WRTETP to write [F^ ] ,p and 
[F2] IP into S3L217. 


Figure 15, — Macro Flow Chart of Overlay (L217 3, 4), FLXFRC 




2 


< NPARTS 




{ Rewind and position SDDATA 
before A 1C. 


I Read row I of [AlC] . 


Call WRTETP to write column I 
of [FGUST] into S3L217. 


Figure 15. — (Concluded) 




Table 10. — Routines Called by FLXFRC 


OVEFLAY (L217, 3,4) 
PFOGRAM FLXFRC 


CALCFF VIP+ 

CALCFG 

DELETF+ 

FETAD+ 

FFTDEL+ 

FNDAIC 

FSF+ 

FSF+ 

INITIR+ 

PFINTM+ 

PEADTP+ 

STAFTR+ 

WFTETP+ 


♦ indicates a routine from the DYLOFLEX 


alternate subroutine library 


) 



Figure 16. — Macro Flow Chart of Overlay (L217, 3, 5), FLXODR 
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Figure 16. — (Continued) 







Table 11. ^ Routines Called by FLXODR 


OVEFLAY(L217,3,5) 
PPOGRAt4 FLXODR 


BINSPT 

DELETR+ 

FETAD+ 

FETDEL+ 

INITIR+ 

PFINTM + 

PFNTCM+ 

READTP+ 

SELPOW 

SFTMFL 

SLCTEL 

STAPTR + 

WRTETP+ 


+ indicates a routine in the DYLOFLEX 
alternate subroutine library 



S3L217 The matrices 

[Ft] (partitioned if necessary) 

[F 2 ] (partitioned if necessary) 
[FG] (by columns) 

FLXODR writes onto the file: 


EOMLOD The following matrices for each surface in a separate logical file 
[X.Y] 1 


{area} i section 6 . 5.2 

[Fi] r 

of volume I 

[P 2 ] 

[FG] 

3.3.6 OVEBLAY(L217,3,6), FLXEOM 

Purpose of FLXEOM 

iTLXEOM generates the file EOMTAP containing the equations of motion matrices 
Ml], [M 2 ], [Mg], [M 4 ], [M 5 ], and [?]. The structural equations of motion ([Mj], 
M 2 ], and [MgJ) are copied from S2L217 where they were stored by STRUCT. TTie 
aerodynamic equations of motion are formed 

IM 4 ] = -[0Fp]'F[SF][F,J 
[M5I = -[^Fpl'T(SFi[F2l 
1 = [^>Fpl%F][FGl 

using the [Fj], [F 2 ], and[PG] stored on S3L217 by FLXFRC. 

The macro flow chart of this overlay is shown in figure 17. 

The routines called are displayed in table 12. 

I/O Devices of FLXEOM 


FLXEOM reads from the files: 


S1L217 The matrices 

(*] 

[<^FP] 


from the first logical file 
from the second logical file 


S2L217 The matrices 

[Ml] 

[M 2 ] (optional) 

[M3] 
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ENTRY 



Figure 17. — Macro Flow Chart of Overlay (L217, 3, 6), FLXEOM 








Table 12. — Routines Called by FLXEOM 


OVEFLAY (L217, 3,6) 
PROGRAM FLXEOM 


CPYM13 


< FETAD+ 
FETDEL+ 
PRINTM+ 
READTP+ 
WFTETP+ 


DELETE + 

FETAD+ 

FETDEL+ 

GENEOM VIP+ 
INITIR+ 

PRINTM+ 

PRNTCM+ 

READTP+ 

SCLMLT 

SETMFL 

STARTR+ 

WRTETP+ 


+ indicates a routine in the DYLOFLEX 
alternate subroutine library 



S3L217 The matrices 


[Fi] 

(partitioned if necessary) 

[Fz] 

(partitioned if necessary) 

[FG] 

(by columns) 


See section 3.6.3 for a detailed description of S1L217, S2L217, and S3L217. 


FLXEOM writes the following matrices onto the file EOMTAP: 


{HEADER} 

[Ml] 

[Mg] 

[M 3 ] 

{FREQM} 

[M 4 ] 

[Ms] 

ik) 

[?] 


(optional) 

See section 6.5.2 
of volume I 


3.4 OVERLAY(L217.4,0), DLAIC 


Purpose of DLAIC 

DLAIC generates aerodynamic force matrices and equations of motion coefficient 
matrices based upon the L216 (DUBFLX)(ref. 4) Doublet Lattice AIC matrices. 


Analytical Steps of DLAIC 
DLAIC performs its task in the following steps: 


1. Print the message NOW IN DLAIC TO..., initialize the count of errors, and 
initialize the scratch file named S1L217. 


2. Call overlay (L217,4,l), DLACDS, to read card set 3.0. 

3. Call overlay (L217,4,2), DLAGEO, to read the Doublet Lattice geometry data from 
NGETP, combine it with card input data from DLACDS, and prepare geometry 
arrays needed in steps 4 through 10. 

4. Call overlay (L217,4,3), DLAMOD, to interpolate for mode shapes at force and 
control points. 


Repeat steps 5 through 8 for each frequency at which the aerodynamics are 
defined. 
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5. Call overlay (L217,4,4), CPANDW, to generate the arrays [CPsb]i [CPGsb]. 

[Wxb]. and [WGtb]- 

6. Call overlay (L217,4,5), FORMWS, to generate the arrays [W] and [WG]. 

7. Call overlay (L217,4,6), SOLVDW, to solve the equation for [Acp], 

[D][Acp] = [WAVG]. 

8. Call overlay (L217,4,7). DLAFRC, to generate the arrays [FR] and [FG]. 

9. Call overlay (L217 ,4,108), DLAODR, to sort the aerodynamic forces per surface and 
write them onto the file EOMLOD. 

10. Call overlay (L217,4,llg), DLAEOM, to prepare [M4], [M5], and [J] and generate 
the file EOMTAP. 

11. Return to the calling program. 

The macro flow chart of this overlay is shown in figure 18. 


The routines called are displayed in table 13. 

The equations solved by DLAIC are summarized as follows: 


[0] 


[FR] = [FAC] 


Q[Asb] [CPsbI 


+ q[l][ATB/IBHD]-l[WAUGJ 


where: “ “ 




[ 0 ] 


[ 0 ] 


(FHCPsbI 


and; 


fCPsBl~2JT [R] +^[R] [^CP^SB v'^f'^CplSBj'*' jj'^CplSB ■‘‘^^'^CplSB 


[FG] = [FAC] 


[ 0 ] 


q[AsB] fCPGgg] 


★ O 


-1 


+ q[-2] [Atb/ib] ID]-* (WG^ugI 


where: [WG^uq] = 
and: [CPGgBl =-^ 


[RG] 


TB 


[R] [GCor] 


[ 0 ] 

SB 


[GC] 


TB 


[ 0 ] 


+ [F] [CPGsg] 


I 
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Loop on 
frequencies. 


Figure 18. — Macro Flow Chart of Overlay (L217, 4, 0), DLAIC 





Figure 18. ~ (Concluded) 
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Table 13. - Routines Called by DLAIC 


OVEFLAY (L217, 4,0) 
PFOGFAM CLAIC 


CPANDW 

- 

OVER LAY (L2 17, 4, 4) 

DLACDS 

- 

OVERLAY (L217,4,1) 

DIAEOM 

- 

OVEFLAY (L217,4,11) 

DLAPRC 

- 

OVER LAY (L2 17, 4, 7) 

DLAGEO 

- 

OVERLAY (L217,4,2) 

DLAMOD 


OVERLAY (L217,4,3) 

DLAODR 

- 

OVERLAY (L217,4,10) 

FORMWS 

- 

OVERLAY (L217,4,5) 

SOLVDW 

- 

OVERLAY (L217,4,6) 

FETAD+ 



FETDEL* 




♦ indicates routines from the DYLOFLEX 


alternate subroutine library 



= '™iReal 
^^2^ ~ Imaginary 

IQR] = (0fp1to/sb I®*'' l'’Rl 
[M 4 ) =-[QRJr,„ 

[M5] Imaginary 

[0] = [0Fp] to/SB 

I/O Devices of DLAIC 

DLAIC prints two messages on the output file; NOW IN DLAIC..., and DLAIC IS 
FINISHED... . All other file communication in DLAIC is performed by its secondary 
overlays. Figure 18 shows the files providing input and output to DLAIC’s overlays. 

In subsections 3.4.2 through 3.4.9 the matrices read from and written onto each file are 
listed. (See secs. 6.4 through 6.5 of vol. I, Engineering and Usage^ for detailed 
descriptions of the format and content of the files for external communication.) Section 
3.6.3 of this volume contains the same infomation for L217 (EOM)’s temporary 
scratch files. 


For a quasi-steady 
solution, the two 
expressions of 

-^in [CPgg] become 
simply i 


3.4.1 OVERLAY(L217,4,l), DLACDS 

Purpose of DLACDS 

DLACDS reads card input data directing the processing of Doublet Lattice (L216) AIC 
matrices. (See card set 3.0 in vol. I of this document, Engineering and Usage,) Options, 
constants, and file names are stored in labelled common blocks for later use. The card 
data defining the correspondence between Doublet Lattice bodies (thin bodies and 
slender bodies) and the L215 (INTERP) modal interpolation surfaces are stored in the 
arrays {IBDTAB} and {ICSTAB} and written onto the temporary storage file S1L217. 

The macro flow chart of this overlay is shown in figure 19. 

The routines called are displayed in table 14. 

I/O Devices of DLACDS 

DLACDS reads card set 3 from the input file and prints the card images on the output 
file with descriptive comments and any necessary diagnostics. 

DLACDS writes on the file: 


S1L217 


The arrays {IBDTAB} and {ICSTAB} (see sec. 3,6.3 for contents) 











Table 14. — Routines Called by DLACDS 


OVEFLAY (L217, 4,1) 
PFOGFAM DLACDS 


FNDKEY+ 

IRDCRD+ 
IFQL + 

NAMFIL* 

WRTETP+ 


+ indicates a routine from the DYLOFLEX alternate subroutine library 
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3.4.2 OVERLAY(L217,4,2), DLAGEO 


Purpose of DLAGEO 

DLAGEO reads the Doublet Lattice (L216) geometry data from the file NGETP, 
combines it with the card input data read from S1L217, and prepares geometry arrays 
governing the calculation of aerodynamic forces. The arrays prepared are written onto 
the first logical file of S1L217 (see sec. 3.6.3). 

Analytical Steps of DLAGEO 
DLAGEO performs its tasks in the following steps: 


1. Initialize geometry arrays and counters. 

2. Call RDLGEO to read the Doublet Lattice (L217) geometry data from the file 
NGETP. 


Rewind S1L217, then repeat steps 3 through 5 for each body (thin or slender) 
defined via card input to DLACDS; IBOD = 1, NBOD. 

3. Read from S1L217 the array {IBDTAB}ibod» card input data associated with 
the body. If control surfaces are attached to the body, the array {ICSTAB} is read 
from S1L217 to define the aerodynamic elements on each of the control sufaces. 

4. Add to the array {ISURF} and build the IBOD portion of the array {IPNSRF}. 
{ISURF} a table of modal interpolation surfaces 

{IPNSRF} a table providing a correspondence between aerodynamic elements 
and modal interpolation surfaces 


5. Build the IBOD portion of the arrays: 

o 

[FAC] Used to double the force of slender body elements on the plane of 

symmetry 

o 

[SF] Used to double the force on elements off the plane of symmetry 

{IGUSTZ} Indicates the gust zone in which each aerodynamic element falls 
[RG] The gust rotation matrix 


6. Build the array [SIGMA] to define the correspondence between interference body 
elements and slender body segments. The [SIGMA] is used to add the interference 
body forces to the slender body s*egments. 

7. Call PRTGFP and PRTGCP to print Doublet Lattice (L216) force point and control 
point geometry. 

8. Transform the reference axis coordinates of all bodies to local structural 
coordinates using the subroutine GTOLT. 
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9. Compress the geometry arrays from NTOTF to NAFP elements by eliminating all 
interference body element entries. 


NTOTF = Total number of force points 

NAFP = Number of aerodynamics force points on all slender and thin 

bodies 

10. Call WRTGEO to write geometry arrays onto the file S1L217 (see sec. 3.6.3). 

11. Read the aerodynamic control matrix from NAETP and extract the reduced 
frequencies. 

12. Return to the calling program. 

The macro flow chart of this overlay is shown in figure 20. 

The routines called are displayed in table 15. 

I/O Devices of DLAGEO 

DLAGEO reads from the files: 

NGETP Doublet Lattice (L216) geometry arrays (see sec. 6.4.4 of vol. I) 

NAETP Doublet Lattice (L216) aerodynamic control matrix (see sec. 6.4.5 of 

vol. I) 

S1L217 Arrays {IBDTAB} and {ICSTAB} 

DLAGEO writes onto the files: 

OUTPUT Geometry tables for force and control points 

S1L217 All geometry arrays required by subsequent DLAIC secondary overlays 

(see sec. 3.6.3) 

3.4.3 OVERLAY(L217,4,3), DLAMOD 

Purpose of DLAMOD 

DLAMOD interpolates over the mode shapes defined in an SA array from L217 
(INTERP) to find: 

Modal deflections at force points 

Modal deflections at control points 

Modal slopes at control points 

Second derivatives of modes at the control points 


[<^FP] 

[<#»cp] 

[‘^►'CP] 

[0"CP] 
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Figure 20. — (Concluded) 




Table 15. — Routines Called by DLAGEO 


OVERLAY (L217,4,2) 
PROGRAM DLAGEO 


GTOLT+ 
IRQL + 
PFTGCP 
PRTGFP 
RDLGEO 

WRTGEO 


FETAD+ 

FETDEL+ 

READTP+ 

WRTETP+ 


+ indicates a routine from the DYLOFLEX 


alternate subroutine library 



Analytical Steps of DLAMOD 

DLAMOD calls AINTL, a DYLOFLEX library routine, to perform the interpolation. If 
more than 20 modes are defined, the arrays are generated and written onto the 
temporary file S1L217 in partitions with up to 20 modes per partition. See section 3.6.3 
for the format of S1L217. 

The macro flow chart of this overlay is shown in figure 21. 


The routines called are displayed in table 16. 

I/O Devices of DLAMOD 

DLAMOD reads from the files: 


SATAP 

S1L217 


{SA} arrays defining modal data 
Geometry arrays 


{SURFTB} 
{IPNSRP} ► 
{X} 

{Y} 

{IPNSCP} 

{XCP} 

{YCP} 


See section 3.6.3 


DLAMOD writes onto the file: 


S1L217 Interpolated modal data (see sec. 3.6.3) 

3.4.4 OVERLAY(L217,4,4), CPANDW 

Purpose of CPANDW 
For a particular frequency CPANDW forms the arrays: 

[CPsbI - -2. [(ri +-^IRl][w'CTlsB^V>*CTlsB]''lft{'»CPlsB + ‘!rt^CP^B^ 

For a quasi-steady solution, the factors multiplying the imaginary parts change to 0.0, 
1/v, and 1/v, respectively. 


[CPGsbI =-^[R1 [RGlsBfGCJsB 
[ WtrI * TB f T b] 
[WGtbI =:J-[RGljBf^CJTB 
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Figure 21. — Macro Flow Chart of Overlay (L217, 4, 3), DLAMOD 
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Figure 21. —(Concluded) 




Table 16. — Routines Called by DLAMOD 


OVEFLAY (L217,4,3) 
PPOGRAM DLAMOD 


AINTL+ 

BINSRT 

DEFIVS 

PETAD+ 

FETDEL+ 

INITIR+ 

PPINTM+ 

READTP+ 

SETMFL 

STAPTF+ 

WPTETP+ 


4- indicates a routine from the DYLOFLEX 
alternate subroutine library 



The macro flow chart of this overlay is shown in figure 22. 

The routines called are displayed in table 17. 

I/O Devices of CPANDW 

CPANDW reads from the file: 

S1L217 Geometry arrays 

[GC] 

[bTg] 

[R] 

[R'J 

and the modal data arrays 

[‘#>CP]sB’ ['^CPJsB’ [‘^CpJsB’ I 

J- See section 3.6.3 

[^CpJtB’ WpJxb J 

CPANDW writes onto the file: 

S1L217 The arrays 

[CPsb] 

[CPGsb] 

[Wtb] 

[WGtb] 

3.4.6 OVERLAY(L217,4,5), FORMWS 

Purpose for FORMWS 

FORMWS forms the matrices: 

[Wtr] 

[W] (F][CP^r1 

( 0 ] 

[WGtbI 

[WGl + [FI [CPGsb] 

and writes them onto the file NTP4 to be used as the right hand sides of the ACp 
solution in SOLVDW. 

Hie macro flow chart of this overlay is shown in figure 23. 
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See section 3.6.3 


► See section 3.6.3 

j 


The routines called are displayed in table 18. 




Figure 22. — Macro Flow Chart of Overlay (L21 7, 4, 4), CPANDW 





Figure 22. - (Concluded) 




Table 17. - Routines Called by CPANDW 


OVEFLAY (L217,4,4) 
PPCX3EAM CPANDW 


CALCP 

CALCPG 

CALCW 

CALCWG 

DELETE + 

INITIR+ 

PRINTM+ 

PRNTCM+ 

READTP+ 

SETMFL 

STARTR+ 

WPTETP+ 


+ indicates a routine from the DYLOFLEX 
alternate subroutine library 



ENTRY 
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( Rewind and call READTP to 
read and [CPGgg] * 


Call READTP to read row I 
of IF], 


Call READTP to read row I 
of [W-tbI. 


jCall WRTETP to write (WJ 


4, 5), FORMWS 





Table 18. — Routines Called by FORMWS 

OVEFLAY (L217,4,5) 

PF06FAM FOFMWS 


ACDWTB 

DELETR+ 

FETAD* 

FETDEL+ 

INITIR+ 

MLTSUB 

PRNTCM+ 

READTP+ 

RETURNF* 

SETMFL 

STARTR+ 

WRTETP+ 


* indicates a routine from the DYLOFLEX 


alternate subroutine library 



I/O Devices of FORMWS 

FORMWS reads from the files: 

NAETP Matrix [F] by rows 

S1L217 Matrices [CPsb] and [CPGsb] 

NTP3 Matrices [W] and [WG] after having written by rows 

FORMWS writes onto the files: 

NTP3 Matrices [W] and [WG] by rows 

NTP4 Matrices [W] and [WG] by columns 

3.4.6 OVERLAY(L217,4,6), SOLVDW 

Purpose of SOLVDW 

SOLVDW solves the following equation for [ACp]: 

[D] [ACp] = [W/WG] 

Analytical Steps of SOLVDW 
SOLVDW performs its task in the following steps: 

1. Determine available length of blank common. 

2. Initialize the files NTP3, NTP4, and NTP5. 

3. Call COPYDQ to copy the quasi-inverse of [D] from NAETP to NTP5. 

4. Call FUTSOL to solve the equation. 

5. Close scratch files and return to the calling program. 

The macro flow chart of this overlay is shown in figure 24. 

The routines called are displayed in table 19. 

I/O of SOLVDW 

SOLVDW reads from the files: 

NAETP Quasi-inverse of [D] 

NTP4 Right-hand-side matrices [W] and [WG] by columns 

SOLVDW writes onto the file: 


NTP3 


Solution columns 



END 


Figure 24. — Macro Flow Chart of Overlay (L217, 4, 6), SOL VDW 
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Table 19. — Routines Called by SOL VDW 

OVERLAY (L217,4,6) 

PROGRAM SOLVDW 


COPYDQ 

FETAD+ 

FETDEL+ 

FUTSOL 


REQFL+ 

RETURNF+ 


FETAD+ 

FETDEL+ 

RETURNF* 


* indicates a routine from the DYLOFLEX 


alternate subroutine library 



SOLVDW also writes/reads data on/from the files MT, NIN, NOUT, and NATAPE 
during the solution inside the subroutine FUTSOL. 


MT = 4LFUT1 

NIN = 4LFUT2 

NOUT = 4LFUT3 

NATAPE = 4LFUT4 

3.4.7 OVERLAY(L217,4,7), DLAFRC 


Purpose of DLAFRC 

DLAFRC generates aerodynamic forces based upon Doublet Lattice (L216) AIC 
matrices: 


[FR] = [FAC] 


[FG] = [FAC] 


[ 0 ] 

+ <j[S ] [Ajg^jg] [ACP] 

q[AsBl [CPsb] 

[ 0 ] ^ = 

+ qf^l tAjByjg] [ACPq] 

icpcgfi] 


Analytical Steps of DLAFRC 

DLAFRC performs its task in the following steps: 

O 

1. Initialize VARDIM storage and read {A}, [£], [FAC], and [CPgs] into core. 


Repeat step 2 for each mode. 

2. Read a column of [ACp] from NTP3, call CALCFS to form a column of [FR], and 
write partitions of [FR] (with IPSIZE columns) onto S3L217. 

3. Initialize {FG}ig and read [CPGsb] from S1L217. 

Repeat step 4 for each gust zone. 

4. Read a column of [ACpc] from NTP3, call CALCFS to form a column of [FG], and 
write (FG}zone onto S3L217. 

5. Return to the calling program. 


The macro fiow chart of this overlay is shown in figure 26. 


The routines called are displayed in table 20. 
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Figure 25. — Macro Flow Chart of Overlay (L217, 4, 7), DLAFRC 
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Table 20. — Routines Called by DLAFRC 


OVEFLAY (L217,4,7) 
PPOGRAM DLAFRC 


CALCFS 

DELETE + 

PETAD+ 

FETDEL+ 

PPNTCM+ 

INITIR+ 

READTP+ 

RETURNF+ 

SETMFL 

STARTR+ 

WRTETP+ 


+ indicates a routine frpm the DYLOFLEX 


alternate subroutine library 


I/O Devices of DLAFRC 


DLAFRC reads from the files: 

S1L217 Matrices {A}, [2], [FAC], [CPsb]> [CPGsb] 

NTP3 Solution columns, [ACp] and [ACpoj 

DLAFRC writes onto the Hie: 

S3L217 Matrices 

[FR] in partitions of IPSIZE columns 

[FG] as NGCP column matrices 
3.4.8 OVERLAV(L217,4,10), DLAODR 

Purpose of DLAODR 

DLAODR sorts the aerodynamic force matrices according to the interpolation surface 
numbers specified for each element via card input data. The resulting surface force 
matrices are written onto EOMLOD. 

The macro flow chart of this overlay is shown in figure 26. 

The routines called are displayed in table 21. 

I/O Devices of DLAODR 

DLAODR reads from the files: 

S1L217 Matrices 

{SURFTB} 

{IPNSRF} 

{X} 

{Y} 

{A} 

S3L217 Matrices [FR] and [FG] (see sec. 3.6.3) 

DLAODR writes on the file: 

EOMLOD For each surface the matrices 

[X.Y]is 

{A}is See section 6.5.2 

[Fi] of volume I 

[F2] 

[FG] ^ 
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Figure 26. — Macro Flow Chart of Overlay (L217, 4, 10), DLAODR 



2 



IP=1NPARTS 

jcati READTP to read (FR) |p. 


S3L217 


SELCRW selects 
rows of (FB) for 
IF, I and (Fjl 


<NPABTS 



NPARTS 


SCLMLT 
IFj) -I/WIF2I 



EOMLOD {call WRTETP to Mirite |F, I and IF 2 I . 


J- t.NGCP 


S3L217, jcall READTP to read |FG|j. 


SELCRW selects 
rows of {fg} j 
for IFGUSTI 


V <NGCP 

V 1 ^ 


<NFREQM 


0 


0 


< NSURFO 



NGCP 

EOMLOD |caM WRTETP to writef GUS^ 


NFREQM 

End file 
EOMLOD 


^ NSURFO 
C END ) 


Figure 26. —(Concluded) 
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Table 21. — Routines Called by DLAODR 
OVEFLAY(L217,4,10) 

pfog’bam dlaodp 


BINSRT 

DELETF+ 

FETAD+ 

FETDEL+ 

INITIR+ 

PPINTM+ 

PRNTCM+ 

READTP+ 

SCLMLT 

SELCPW 

SELROW 

SETMFL 

SLCTEL 

STARTR+ 

WPTETP+ 


♦ indicates a routine from the DYLOFLEX 


alternate subroutine library 



S.4.9 OVERLAY(L217,441), DLAEOM 

Purpose of DLAEOM 

DLAEOM generates the aerodynamic equations of motion matrices ([M 4 ], [M 5 ], and 
[^]) and Mrrites the file EOMTAP. The matrix equations are: 

[M4] = -[0Fpl^fSF](FR]Rg3j 

IM5I = -[#Fpl’'lSFHFRIlmwnary 

[?] = 1<^FP1^[SF]IFG] 


The macro flow chart of this overlay is shown in figure 27. 
The routines called are displayed in table 22 . 

I/O Devices of DLAEOM 


DLAEOM reads from the files: 

S1L217 Matrices [SF] (from the first logical file) and [d>Fp] (from the second 

logical file) 

S3L217 Matrices [FR] (by partition) and [FG] (by columns) 

DLAEOM writes onto the file EOMTAP the matrices: 

{HEADER} 

[Ml] 

[M 2 ] (only if damping requested) 

[M3] 

{FREQM} 

Repeat [M 4 ] 

per frequency < [M 5 ] 

{fj^} (only for first frequency) 

3.6 OVERLAY(L217,6,0), DLPRES 
Purpose of DLPRES 

DLPRES prepares the equations of motion file, EOMTAP, and the aerod 3 rnamic force 
file, EOMLOD (sec. 3.6.2). It uses aerodynamic information from Doublet Lattice (L216) 
(ref. 4) the generalized forces and the unsteady aerodynamic pressure. 
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Figure 27. —Macro Flow Chart of Overlay (L217,4, 11), DLAEOM 





Call READTP to read {FG} 


Call WRTETP to write 
[M4] and [Mg] . 


Figure 27. — (Continued} 






Table 22. — Routines Called by DLAEOM 


OVERLAY (L2 17, 4, 11) 
PFOGPAK DLAEOM 


CGNEOM 


CYPMI 3 


FETAD+ 

• FETDEL+ 
PRINTM+ 
READTP+ 
^ WRTETP+ 


DELETP+ 

FETAD+ 

FETDEL+ 

INITIR+ 

PRINTM+ 

PRNTCM+ 

READTP+ 

SCLMLT 

SETMPL 

STARTR+ 

WRTETP+ 


♦ indicates a routine from the DYLOFLBX 
alternate subroutine library 



Analytical Steps of DLPRES 

DLPRES performs its task in the following steps: 

1. Print the message NOW IN DLPRES TO..., initialize the count of errors, and 
initialize the scratch file named S1L217. 

2. Call overlay (L217,5,l), DLPCDS, to read the card input data directing the 
processing of Doublet Lattice (L216) generalized forces and unsteady pressures. 

3. Call overlay (L217,5,2), DLPGEO, to process the Doublet Lattice (L216) geometry 
data and combine it with the card input date to prepare arrays directing the 
processing of steps 4 through 6. 

4. Call overlay (L217,5,4) to generate the aerodynamic forces per reduced frequency 
and per mode. 

6. Call overlay (L217,5,5) to sort the aerodynamic forces and write them onto 
EOMLOD by surface, 

6. Call overlay (L217,5,6) to generate the equations of motion matrices and write 
them onto EOMTAP. 

7. Close the scratch file S1L217 and return to the calling program. 

The macro flow chart of this overlay is shown in figure 28. 

The routines called are displayed in table 23. 

The equations solved by DLPRES are summarized as follows: 


IFR] = IFAC][S]]A][ACP] 


IFi] = q[FRjR^^ 


_q_ 

(^ 2 ^ “ CO f^^hmaginary 

(ifco = 0 , [F2] =[ 0 j) 

]FG! =^{FRln 

(if <0 - 0 , -L rep laces -I3.) 
V cu 


where n indicates column n of [FR] 

n = NTDEF ifcJ>0 
n=NPDEF ifcj = 0 


[M4] --2q[QlReal 

[M 5 ] =-^tQ]imaginary (if = 0, (M 5 ] = 10] ) 

{ ^ ^ {Qlmaginary(NTDEF " ^ {^Reall NTDEF 

{^} NPDEf} (when« = 0) 
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Figure 28. —Macro Flow Chart of Overby (L217, 5, 0), DLPRES 




Table 23. — Routines Called by DLPRES 
OVERLAY (L217;,5,0) 


PROGRAM 

DLPRES 

DLPCDS 

OVERLAY (L2 17, 5,1) 

DLFEOM 

OVERLAY (L217,5,6) 

DIPFRC 

OVERLAY (L2 17, 5, 4) 

DLPGEO 

OVERLAY (L21 7, 5, 2) 

DLPODR 

OVERLAY (L217,5,5) 

FETAD+ 


FETDEL+ 



+ indicates routines from the DYLOFLEX 


alternate subroutine library 


I/O Devices of DLPRES 


DLPRES prints two messages on the output file; NOW IN DLPRES..., and DLPRES IS 
FINISHED... . All other file communication in DLPRES is performed by its secondary 
overlays. Figure 29 shows the files used in DLPRES' secondary overlays. 

In subsections 3.5.1 through 3.5.5 the matrices read from and written onto each file are 
listed, (See sec. 6.4 to 6.5 of vol. I, Engineering and Usage, for detailed descriptions of 
the format and content of files for external communication.) Section 3.6.3 of this volume 
contains the same information for L217 (EOM)'s temporary scratch files. 

3.5,1 OVERLAY(L217,5,l), DLPCDS 

Purpose of DLPCDS 

DLPCDS reads card input data directing the processing of Doublet Lattice (L2I6) 
generalized forces and unsteady pressures. (See card set 3.0 in vol. I of this document, 
Engineering and Usage.) 


Analytical Steps of DLPCDS 
DLPCDS performs its task in the following steps: 

1. Initialize options and constants to their default values. 

2. Read the constant card and extract RHO, VEL, NTDEF, and NPDEF. 

Repeat step 3 until an END card is read. Then, jump to step 4. 

3. Read another data card, extract the keyword, and jump to the location which 
processes that card. The keywords are: 


SATAP 

CHECK 

SLENDER 

NABTP 

PRINT 

THIN 

NGETP 

SELECT 

CONTROL 

EOMTAP 


END 

EOMLOD 




• The keyword card CONTROL is recognized only after the THIN card. 

• The information from the THIN or SLENDER card is stored in the array 
IBDTAB and written onto the file S1L2I7. 

• The panels specified after a CONTROL card are stored in the array ICSTAB 
and written onto the file S1L217. 

4. Write and end of file onto the file S1L217. 
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Figure 29. —Files Providing Input/Output to DLPRES Overlays 




5. Return to the calling program. 


The macro flow chart of this overlay is shown in figure 30. 

The routines called are displayed in table 24. 

I/O Devices of DLPCDS 

DLPCDS reads all of the card input data directing the processing of Doublet Lattice 
(L216) generalized forces and unsteady pressures. 

It prints all of the card images, additional descriptive comments, and diagnostics for all 
errors discovered. 

DLPCDS writes arrays of surface data onto the file S1L217. The arrays are {IBDTAB} 
and {ICSTAB}. See section 3.6.3 for the exact format and content of S1L217. 

3.5.2 OVERLAY(L217,5,2), DLPGEO 

Purpose of DLPGEO 

DLPGEO reads the Doublet Lattice (L216) geometry data from the file NGETP, 
combines it with the card input data read from S1L217, and prepares geometry arrays 
governing the calculation of aerodynamic forces. The arrays prepared are written onto 
the first logical file of S1L217 (see sec. 3.6.3). 

Analytical Steps of DLPGEO 

DLPGEO performs its tasks in the following steps: 

1. Initialize geometry arrays and counters. 

2. Call RDLGEO to read the Doublet Lattice (L216) geometry data from the file 
NGETP. 

Rewind S1L217, then repeat steps 3 through 5 for each body (thin or slender) 
defined via card input to DLPCDS; IBOD = 1, NBOD. 

3. Read from S1L217 the array {IBDTAB}ibqd, the card input data associated with 
the body^ If control surfaces are attached to the body, the array {ICSTAB} is read 
from S1L217 to define the aerodynamic elements on each of the control surfaces. 

4. Add to the array {ISURF} and build the IBOD portion of the array {IPNSRF}. 
{ISURF} a table of modal interpolation surfaces 

(IPNSRF) a table providing a correspondence between aerodynamic elements 
and modal interpolation surfaces 
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Table 24. — Routines Called by DLPCDS 

OVERLAY (L2 17,5, 1) 

PRORAM DLPCDS 


FNDKEY+ 

IRDCRD+ 

NAMFIL+ 

WRTETP* 


♦ indicates a routine from the DYLOFLEX 


alternate subroutine library 


5. Build the IBOD portion of the arrays: 

O 

[FAC] Used to double the force of slender body elements on the plane of 

symmetry 

O 

[SF] Used to double the force on elements off the plane of symmetry 

{IGUSTZ} Indicates that all elements are in gust gradual penetration zone 

number one. Gradual penetraton is not allowed in DLPRES. 

6. Build the array [SIGMA] to define the correspondence between interference body 
elements and slender body segments. The [SIGMA] is used to add the interference 
body forces to the slender body segments. 

7. Call PRTGFP to print the Doublet Lattice (L216) force point geometry. 

8. Transform the reference axis coordinates of all bodies to local structural 
coordinates using the subroutine GTOLT. 

9. Compress the geometry arrays from NTOTF to NAFP elements by eliminating all 
interference body element entries. 

NTOTF = Total number of focal points. 

NAFP = Number of aerodynamic force points on all slender and thin bodies. 

10. Call WRTGEO to write geometry arrays onto the file S1L217 (see sec. 3.6.3). 

11. Return to the calling program. 

The macro flow chart of this overlay is shown in figure 31. 

The routines called are displayed in table 25. 

I/O Devices of DLPGEO 

DLPGEO reads the Doublet Lattice (L216) geometry data from the file NGETP (see 
sec. 6.4.4 of vol. I). 

DLPGEO reads the arrays {IBDTAB} and {ICSTAB} from S1L217. They contain the 
card input data read by DLPCDS. 

DLPGEO prints the geometry arrays and then writes them onto the file S1L217. 

3.5,3 OVERLAY(L217,6,4), DLPFRC 

Purpose of DLPFRC 

DLPFRC calculates the complex aerodynamic forces on the structural elements (the thin 
body boxes and slender body segments). Note that the interference body forces are 
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Figure 31. -^Macro Flow Chart of Overlay (L217, 5, 2), DLPGEO 




m 




< NBOD 



Figure 31. —(Concluded) 











Table 25. — Routines Called by DLPGEO 

OVERLAY (L217, 5,2) 

PROGRAM DLPGEO 


GTOLT* 

IRQL+ 

PRTGFP 

PDLGEO 


WRTGEO 


FETAD* 

FETDEL* 

READTP+ 

WRTETP+ 


♦ indicates a routine from the DYLOFLEX alternate subroutine library 
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included in the slender body segment forces by the summation matrix [S] in the 
following equation: 

[Fi2]j^=[FAC][i][A][ACP],, 

k = 1, NFREQM (the number of reduced frequencies) 

Analytical Steps of DLPFRC 
DLPFRC performs its task in the following steps: 

1. Initialize VARDIM storage, arrays, and counters. 

2. Read [A], [X], and [FAC] from S1L217. 

3. Read the {ACM}, the aerodynamic control matrix, from NAETP and extract the 
reduced frequencies. 

Repeat steps 4 through 7 for each reduced frequency. 

4. Skip the [F] and [D] quasi-inverse matrices if they are on NAETP. 

Repeat steps 5 through 7 for each mode; IDOF = 1, NDOF. 

5. Read {ACP)idof from NAETP. 

6. Call subroutines CEQAC, SUMFRC, and CEQAC again to form 

[A]|ACP|. 

7. Call WRTETP to write (Fi2}iD 0F onto S3L217. 

8. Close the scratch files and return to the calling program. 

The macro flow chart of this overlay is shown in figure 32. 

The routines called are displayed in table 26. 

I/O Devices of DLPFRC 

DLPFRC reads from the files: 

NAETP {ACM} Aerodynamic control matrix 

[ACP] Complex unsteady pressure matrix by columns. 

(See sec. 6.4^5 of vol. I and ref. 4 for a detailed description of NAETP contents.) 
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ENTRY 


initialize VARDIM 
storage, the files 
NAETPand S3L217, 
and all arrays used 


Rewind 

NAETP 

S1L217 

S3L217 


S1L217 


NAETP 


Call READTP to read 

[A], [2], and [FAc]. 


Call READTP to read 


Extract (kVALUES} 
from {ACM} and 
convert to { CPS } 


IFREQ=1,NFREQM jfrepSencies. 


lAERO 


SKIPFD skips 
[F] and [D]q‘ 



IDOF = 1, NDOF |Loop on modes. 


NAETP {Call READTP to read 


)4CP| 


Figure 32. —Macro Flow Chart of Overlay (L217, 5, 4), DLPFRC 





Table 26. — Routines Called by DLPFRC 


OVERLAY (L2 17, 5,4) 
PROGRAM DLPFRC 


CEQAC 

FETAD+ 

FETDEL+ 

INITIR+ 

READTP+ 

SETMFL 

SKIPFD 

STARTR+ 

SOMFRC 

WRTETP+ 


+ indicates a routine from the DYLOFLEX 
alternate subroutine library 



S1L217 

[A] 

Aerod 3 mamic element areas 


[2] 

Force summation matrix 


[FAC] 

Array to double the forces on slender body elements on the plane 
of symmetry 

DLPFRC 

writes onto 

the file: 

S3L217 

[F 12 ] 

Unsteady aerodynamic forces on thin body boxes and slender body 
segments 

Matrix is written by columns and repeated per reduced frequency 


3.5.4 OVERLAY(L217,5,5), DLPODR 

Purpose of DLPODR 

DLPODR sorts the aerodynamic force matrices according to the interpolaton surface 
numbers specified for each element via card input data. The resulting surface force 
matrices are written onto the file EOMLOD. 

The macro flow chart of this overlay is shown in figure 33. 

The routines called are displayed in table 27. 


I/O Devices of DLPODR 

DLPODR reads from the file: 

S1L217 {SURFTB} Interpolation surface table 

{IPNSRF} Panel/surface correspondence table 

{X} X-coordinates of aerodynamic force points 

{ Y} Y-coordinates of aerod 3 mamic force points 

S3L217 {a} Aero element areas 

[F 12 ] Complex aerodynamic forces at thin and slender body force 
points 


DLPODR writes on EOMLOD for each interpolation surface IS the following matrices: 

[X,Y]is Liocal coordinates of the surface’s force points 
{A}is Areas of surface IS elements 

[F|] Real part of aero forces on surface IS 

[F 2 ] Imaginary part of aero forces on surface IS 

[FG] Complex matrix of gust forces on surface IS 
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Figure 33. -Macro Flow Chart of Overlay (L217. 5. 5), DLPODR 







Figure 33. — (Concluded) 





Table 27. — Routines Called by DLPODR 

OVERLAY (L2 17, 5, 5) 

PROGRAM DLPODR 


BINSRT 

DELETE + 

FETAD+ 

FETDEL+ 

INITIR+ 

PRINTM+ 

PRNTCM+ 

READTP* 

SCLFRC 

SELCRW 

SELROW 

SETMFL 

SLCTEI 

STARTR+ 

WRTETP+ 


♦ 


indicates routines from the DYLOFLEX 
alternate subroutine library 



3.6.5 OVERLAY(L217,5,6), DLPEOM 

Purpose of DLPEOM 

DLPEOM generates the file EOMTAP containing the equations of motion coefficient 
matrices. The aerodynamic matrices ([M 4 ], [M 5 ], and [d>]) are based upon the Doublet 
Lattice (L216) (ref. 4) generalized forces. 

[M4] =-2{QReall 
2 

[M5] = ^(QlmaginaryJ 

[0] =-^{QReal}NPDEF<^'^^®" 

{^Imaginary I nTDEF " ' {QReal jNTDEF,) <'^^®" " * 

where [Q] = the generalized air force matrix from Doublet Lattice (L216) 

The macro flow chart of this overlay is shown in figure 34, 

The routines called are displayed in table 28. 


I/O Devices of DLPEOM , 

DLPEOM reads from the files: 

NAETP [Q] Complex matrix of generalized forces from Doublet Lattice 

(ref. 4) 


S2L217 [Ml] 
[Ms] 


Generalized stiffness matrix 
Generalized damping matrix (optional) 
Generalized mass matrix 


DLPEOM writes onto the file 


EOMTAP 


{HEADER} 

[Ml] 

[M2] 


DYLOFLEX header matrix introducing following arrays 
From S2L217 above 


[Ms] 

(FREQM) Frequencies at which aerodynamics are defined 


{M 4 } 

{Mg} 

Aerodynamic force matrices 

{h} 

Distance to gust zones 

[5] 

Gust force matrix 
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tnitifliize VAR DIM 
storage and the 
file EOMTAP 



Figure 34. —Macro Flow Chart of Overlay (L217, S, 6), DLPEOM 




Figure 34. —(Concluded) 
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Table 28. — Routines Called by DLPEOM 


OVERLAY (L2 17, 5, 6) 
PROGRAM DLPEOM 


CPYM13 


. PET AD ♦ 
FETDEL+ 
PRINTM+ 
READTP+ 
WRTETP* 


DELETR+ 

PETAD+ 

FETDEL+ 

INITIR+ 

PRINTM+ 

PRNTCM+ 

READTP+ 

SCLPHT 

SETMFL 

SKIPFD 

SPLITQ 

STARTF+ 

WRTETP+ 


+ Indicates a routine from the DYLOFLEX 
alternate subroutine library 


3.6 DATA BASES 


The program’s data bases include input and output files plus internal scratch 
(temporary) storage files and common blocks. 

8.6.1 INPUT DATA 

Data is input to the program in two forms; data cards and magnetic files. 

Card Input Data 

For a complete description of L217 (EOM) card input data, see section 6.3 in volume I of 
this document. 

Magnetic File Input Data 

For a complete description of L217 (EOM) disk or tape input data, see section 6.4 in 
volume I of this document. 


The files read are: 




KTAPE(GSTIF) 


NGETP 

from L216 (DUBFLX) 

MTAPE(GMASS) 


GDTAPE 

from FLEXSTAB 

SATAP 

from L215 (INTERP) 

SDSSTP 

from FLEXSTAB 

NAETP 

from L216(DUBFLX) 



3.6.2 OUTPUT DATA 





The results of a L217 (EOM) run will be printed, and written on magnetic files. 

Printed Output Data 

For a complete description of the printed output data, see section 6.5.1 in volume I of 
this document. 

Magnetic File Output Data 

For a complete description of L217 (EOM) disk or tape output data, see section 6.5.2 in 
volume I of this document. The files written are: 

EOMTAP 

EOMLOD 

8.6.3 INTERNAL DATA 

Two methods are used to pass data firom one portion of the program to another. They are 
labelled common blocks and scratch (temporary) storage disk files. 
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Magnetic Disk Files 

L217 (EOM> uses 10 disk files for temporary (scratch) storage of data. The flies are: 

S1L217, S2L217 
S3L217 

NTP3, NTP4, NTP5, MT, NIN, NOUT, and NATAPE 

The first two files are used in all L217 primary overlays. The S3L217 is used in all 
aerodynamic primary overlays. The other files are used only in DLAIC when solving the 
system of complex simultaneous equations. All scratch files are released before L217 
(EOM) terminates execution. 

The first three files are always written with the subroutine WRTETP and read with the 
subroutine ^ADTP. The READTPAVRTETP format is described in reference 1. 1%e 
other files are written with unformatted FORTRAN write statements. 

S1L217 

The scratch file named S1L217 may contain several logical files of data, each 
terminated by. an end of file (EOF). The data on S1L217 is always in the 
READTP/WRTETP format. The matrices existing on the file vary according to which 
overlay was executed last (see figs. 35 through 39). 

Ihe following matrices are written onto S1L217 inside STRUCT’s secondary overlays: 

[CSDATA] Data for a control surface adding freedoms. Each row of CSDATA 

represents a control suface node and contains six numbers as input 
on card 2.9.1. . : 


Y f Coordinates of control surface node 


. m = Mass 

ly = Moment of inertia about local y*axis 
Ix - Moment of inertia about local X'axis 

[A] =i^cs1^imcs][^st1 

[C] =[0csl^lMCS][^Cs^ 

All of the data listed in figure 35 is released just before overlay (L217,2,0), STRUCT 
terminates execution. 

The S1L217 written by overlays (L217,3,l), FLXCDS 

(L217,4,l), DLACDS 
(L217,5,l), DLPCDS 
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Matrix Size Matrix 



Original generalized mass matrix 
written in RSTRCT only if 
reat^ from cards « 

Original generalized stiffness 
matrix written in RSTRCT only 
if read from cards. 


Damping factors written in 
RSTRCT only if read from cards. 


Written in RSTRCT only if cards 
specifying control surface 
freedoms were read. 


NCSEF = number of control 
surfaces adding 
freedoms 


Written in FORMAC only if 
control surface freedoms 
added . 


Figure 35. — S1L217 at the End of Overlay (L21 7,2,3,), Form 13 
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Ma t rix 


NSl®F X 
NAFP X 
NAFP X 
NAFP X 
NACP X 
NACP X 
NACP X 
NTOTF X 
NTOTF X 
NAFP X 

NACP X 

NAFP X 

NACP X 
NSBFP X 
NSBFP X 


D escription 

Surface Table 
Force Point/Surface Table 
Force Point x-coordinates 
Force Point y- coordinates 
Control Point/Surface Table 
Control Point x-coordinates 
Control Point y-coordinates 
Areas of all aero elements 
Force summation matrix ([Z]) 

Array of 1*s and 2’s to double force 
of slender body elements on the P.O.S. 
Array of integer numbers indicating the 
gust zone of each control point. 

Array of 1's and 2's to double force 
of elements off the P.O.S. 

Gust rotation matrix 

Slender body segment midpoint radii 

Slender body segment midpoint radius 

slopes. 


♦ Meaningless with FLEXSTAB data 
♦♦Meaningless except with Doublet 
Lattice (L216) AIC data 

Figure 36. — S1L217 Logical File 1 After a Geometry Overlay 
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M atr i 

NACP 

NACP 

NAFP 


Siz e Ma tri x 


Des cr iption 



Modal deflections at control points 
Modal slopes at control points 
Modal deflections at force points : 

The three matrices are repeated 
per partition of the total 
number of modes with 


IPSIZE < 20 


Figure 37. - S1L217 Logical File 2 After Overlay (L21 7,3,3}, FLXMOD 



£liix:is_si2£ 


Matrix 


De scription 


NSBCP 

NSBCP 

NSBCP 

NTBCP 
NTBCP 
MAFP X 



Modal deflections at SB control points 
Modal slopes at SB control points 
2nd deriative of modal deflections 
at SB control points. 

Modal deflections at TB control points 
Modal slopes at TB control points 
Modal deflections at all force points 

The matrices are repeated 
per partition of the total 
number of modes with 
IPSIZE < 20 


Figure 38. — S1L217 Logical File 2 After Overlay (L217,4,3), DLAMOD 



tiaiL£i2{.§is£ 


De scription 


NSBFP2 X NSBFF 
NSBFP2 X NGCP 
2 X NCOF 
2 X NGCP 
2 X NDOF 
2 X NGCP 


L«tbJ 
Kb J 
KbJ 
LwGtb J 


Slender body pressure coefficients 
Slender body gust pressure coefficients 
Norinalwash on thin body element 1 
Gust normalwash on thin body element 1 

2 

2 


2 X NDOF 
2 X NGCP 


^'^TBjljBpp Normalwash on thin body element NTBFP 
KbLf, Gust normalwash on thin body element NTBFP 


EOF 


Note: 



Because all of these 
arrays are complex the 
number of rows written 
onto S1L217 were doubled. 


F/gure 39. —S1L217 Logical File 3 After Overlay (L217.4,4), CPANDW 



The card input data for all three aerodynamic paths is read from cards and stored on 
S1L217 in the same fashion. The file will contain the array {IBDTAB} for each thin or 
slender body specified by a card. The array {IBDTAB} contains seven elements: 


Description 

Interpolation surface number for the 
body’s Y-motion (in the left most 30 bits) 

Interpolation surface number for the body’s 
Z-motion (in the right most 30 bits) 

Surface number in which the body’s forces 
will be placed when written on EOMLOD 


Geometric shift values 


Number of control surfaces attached to 
the surface (= 0 for slender bodies and 
^10 for thin bodies) 

Body type code 

(1 = slender body and 3 = thin body) 

If NCS > 0, the body table array will be followed immediately by NCS arrays 
(ICSTAB). The {ICSTAB} will contain the numbers of the thin body elements 
contained in the control surface. Note that the elements are numbered relative to each 
thin body, not for the entire structure. 

The S1L217 logical file 1 written by overlays (L217,3,2), FLXGEO 

(L217,4,2), DLAGEO 
(L217,5.2), DLPGEO 

Although the geometry input data for the three aerodynamic paths is different, the 
geometry stored on the first logical file of S1L217 by the geometry processors is in 
exactly the same format. 

Figure 36 displays the matrices written on the file. The following descriptions 
supplement the short ones in the figure. 

{SURFTB} 

Each element of SURFTB is an integer representing a modal interpolation surface used 
in L217 (EOM). 


Element 

1 ISY 

ISZ 

2 ISTORE 

3 XSHIFT 

4 YSHIFT 
6 zsraFT 

6 NCS 

7 IBDTYP 
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SURFTBi « 10® 

•IS 

Interpolation surface number 

+ 10® 

•ISTORE 

Output surface number 

+ 10® 

• ITYPE 

Body type code (1 = SB, 3 = TB) 

+ 10 

•NMOT 

Number of motion types for 
Body (1 = TZ, 2 = TY, 3 = both) 

+ 

IPOS 

Plane of symmetry code 
(0 = off, 1 = on) 


{IPNSRF} and {IPNSCP} 

Each element of IPNSRP/IPNSCP is an integer representing a force/control point. 


IPNSRF 

= 10® 

•IS 


+ 10® 

• ISTORE 


+ 10® 

•ITYPE 


+ 

IPOS 


{RG} 

The gust rotation matrix has elements for all thin body and slender body control points. 

For a thin body element: 

Symmetric analysis RG{ = cosy} 

Anti-symmetric analysis RGi = -sinyi 
Where yi is the element’s dihedral angle 
For a slender body element: 

Symmetric analysis 

Y-motion element RG| = 0.0 
Z-motion element RGi = 1.0 

Anti-symmetric analysis 

Y-motion element RG| = 1.0 
Z-motion element RGi = 0.0 

[SIGMA] 

Each row i of [SIGMA] (also known as [S]> indicates where the forces on aerodynamic 
element i are to be stored. The primary purpose of SIGMA is to add the interference 
body forces to the slender body forces. However, in the following paragraph, other 
factors are included in SIGMA for a FLEX6TAB AIC data case. 
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SlGMA(i,l) - Indicates which force point will receive the SIGMA(i,2) fraction of the force 
on aero element i 

SIGMA(i,3) - Indicates which force point will receive the SIGMA(i,4) fraction of the force 
on aero element i 

For a FLEXSTAB AIC data case, the following ACp factors are included in the second 
and fourth columns of SIGMA: 

7t/2 for slender body Y-motion elements 

"- 7 r /2 for slender body Z-motion elements 

— 1.0 for interference body elements 

1.0 for thin body elements 

The S1L217 logical file 2 written by overlays (L217,3,3), FLXMOD 

(L217,4,3), DLAMOD 

The modal deflections and slopes calculated by FLXMOD and DLAMOD are written 
onto the second logical file of S1L217. 

Overlay (L217,4,4), CPANDW, writes the intermediate arrays onto the third logical file 
ofSlL217 (fig. 39). 

S2L217 

S2L217 is written by overlay (L217,2,0), STRUCT (see fig. 40). It will contain the 
generalized structural matrices after STRUCT is finished. Actually, they are written in 
the secondary overlay, overlay (L217,2,3), FORM13. The matrices are written in the 
READTP/WRTETP format. 

S3L217 

Overlays (L217,3,4), FLXFRC; (L217,4,7), DLAFRC; and (L217,5,4), DLPFRC all write 
the unsteady aerod 3 mamic forces onto S3L217 with the subroutine WRTETP, However, 
the exact content and format is different for each. FLXFRC output is displayed in 
figure 41. In DLAFRC, the entire logical file (including the end of file) is repeated for 
each frequency at which the aerod 3 mamics are defined. Because the force matrices [FR] 
and [FG] are complex, they are written onto S3L217 with NAFP2 rows where 
NAFP2 = 2 * NAFP (see fig. 42). 

In DLPFRC, because the force matrices [Fj2]kJ k = 1, NFREQM are complex, the 
matrix columns are written onto S3L217 with NAFP2 rows where NAFP2 = 2 * NAFP. 
The columns of [Fi2]i are written first followed by the columns of [Fj 2] 2> 6^- See 
figure 43. 

NTP3 

FORMWS calls WRTETP to write the rows of [W], the augmented normalwash, and 
[WG];>the augmented gust normalwash, onto NTP3 (see fig. 44) as they are formulated: 
FORMWS later reads the rows back into core and writes the same data onto NTP4 by 
columns with FORTRAN unformatted write statements. 


5aatrix_Size Matrix 


Descri p tion 


NDOF X NDOF 
1 X NDOF 


NDOF X NDOF 



Generalized stifness matrix 
Diagonal elements of the 
generalized damping matrix. 
Generalized mass matrix 

♦Ma is present if and only if 
structural damping was 
specified with card 2.8. 


Figure 40. -S2L217 After Overlay (L217,2,0). STRUCT 



M atrix Siz e Matrix Descript4.on 


NAFP 

NAFP 

NAFP 

NAFP 


NAFP 



Figure 41. - S3L21 7 After Overlay (L217,3,4). FLXFRC 



Matrix 


Matrix 

NAFP2 

NAFP2 : 
NAFP2 5 

• 

NAFP2 : 


.Size Matrix S®S££iE£i2D 



The aerodynamic forces on 

Thin and Slender Body force points 

Repeat per modal partition; 

IP=1, NPARTS. 

The NGCP columns of 
the gust aerodynamic 
forces on Thin and 
Slender Body force points. 


Figure 42. -S3L217After Overlay (L21 7,4 J), DLAFRC 


£jatr|.jc_Siie Matrix 


NAFP2 X 1 
NAFP2 X 1 


NAFP2 X 
NAFP2 X 
NAFP2 X 


NAFP2 X 


Descripti on 



per generalized 
frequency. 


Figure 43. - S3L217 After Overlay (L21 7,5,4), DLPFRC 



M atrix Size Matrix 


Description 



row 1 
row 1 

2 

2 


row NBOX 
row NBOX 


Figure 44. - NTP3 After Overlay (L217,4,S), FORMWS 
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SOLVDW calls the subroutine FUTSOL to solve the system of equations: 

[D][Acp] =W/WG] 

FUTSOL reads the right hand sides, [WAVG], by columns from the file NTP4. 

FUTSOL reads the D-quasi-inverse matrix in a variable number of records from the 
file NTP5. 

FUTSOL writes the solution, [ACp], by columns onto the file NTP3 (see fig. 45) with 
FORTRAN unformatted write statements. Note that the matrix is complex. 

NTP4 

Overlay . (L217, 4,5) FORMWS uses a FORTRAN unformatted write statement to write 
the columns of [W] and [WG] onto the file named NTP4 (see fig. 46), 

e.g., For column J of [w] it uses: 

WRITE (NTP4) (W(I,J), 1=1, NBOX) 

Note that both matrices are complex. 


NTP5 

Overlay (L217,4,6), SOLVDW, calls the subroutine COPYDQ to read the 
D-quasi-inverse records from the file NAETP, and writes the identical records onto 
NTP5. SOLVDW then calls the subroutine FUTSOL to solve the system of equations: 

[DHACp] = [W/WG] 

Copying the D-quasi-inverse to NTP5 allows FUTSOL to repeatedy reposition the file 
without destroying the sequential processing of the file NAETP. 

For a description of the records comprised in the D-quasi-inverse matrix, see the 
description of NAETP in reference 4. 

MT, NIN, NOUT, and NATAPE 

The scratch files MT, NIN, NOUT, and NATAPE are defined, used, and released within 
the subroutine FUTSOL. No other routines write or read these files. All of them are 
written/read with FORTRAN unformatted write/read statements. 

NIN contains a temporary copy of the solution columns. 

MT, NOUT, and NATAPE contain modified versions of the right hand sides and the D- 
quasi-inverse trapezoidal matrix. 



Matrix Si 2 a Matrix 


Descr iptio n 


NBCX X 1 

1 

{ACP> 1 

AGP due to mode 1 

NBOX X 1 

• 

{ACP>2 

• 

Acp due to mode 2 

• 

• 

• 

NBOX X 1 

• 

• 

Acp due to mode NDOF 

NBOX X 1 
• 

^^^G^NDOF+J 

ACp due to gust on zona 1 
• 

• 

NBOX X 1 

• 

^'^'^G^NDOF 
+NGCP 
BO^, 

• 

• 

ACp due to gust on zone NGCF 


Figure 45. - NTP3 After Overlay (L21 7,4,6), SOL VDW 



Matrix size 

Matrix 

Description 

NBOX X 1 

{W>i 

Augmented normalwash, column 1 

NBOX X 1 
• 

{W>2 

• 

Augmented norma Iwash, column 2 
• 

• 

NBOX X 1 

• 

• 

• 

Augmented norma Iwash, column NDOF. 

NBOX X 1 

• 

{WG>i 

• 

Augmented gust normalwash, column 1. 
• 

• 

NBOX X 1 

• 

• 

• 

Augmented gust normalwash, column NGCP 


EOF 



Figure 46. - NTP4 After Overlay (L21 7.4.5). FORMWS 
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Common Blocks 


Table 29 shows the common blocks used in the program and the overlays in which they 
are defined. 

The labelled common blocks are used for communication between the primary and 
secondary overlays, and for communication between routines in a secondary overlays. 
The block names and contents are described in table 30. 

Blank common is used as a variable length working area. In general, the main program 
of an overlay calculates the area required for arrays in the various subroutines and 
passes a dimension and first word addresses of the arrays through the subroutines’ 
argument lists. 
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Table 29. - Common Blocks Defined in Each Overlay 


COMMON 

BLOCK 


OVERLAY 





CSI 

Pk 

Pk 

Cx 

Pk 


o 

o 

g 

Pk 

P 

< 

A 

A 

A 


S g ^ ! eJ 3 

6 Z Q S CO < Q D 

U O U U O O Q 


CO CO CO 

^ ^ ^ 

!S I SS 

g s ^ 


o CO H G pj M 

H H Z < G hi 

H Z H e H 2 M 

Pk < Z 7 < U H 

O Pk P< Z CO CO CO 



VARDIM 






Table 30. — Contents of Common Blocks 

Each of the common blocks is described on a separate page* The- 
following abbreviations are used to describe the individual 
variables. 

NO. - indicates the variable number in the 

common block list. 

VARIABLE - common block item name. 

T - type of variable 

I = integer 

F = real 

C = complex 

H = hollerith j 

DIMENSION - number of elements in variable. 

ENG. NOM. - engineering nomenclature - symbolic. 
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Table 30, — (Continued) 


LABEi^ED COMMON NAME; BODPAN 


DESCRIPTION: Problem size in number of aero bodies and elements 


NO. 

VARIABLE 

T 

DIMENSION 

ENG. 

NOM. 

DESCRIPTION 

1 

NBOD 

I 



Number of zero bodies (Thin 
and Slender) 

2 

NTB 

I 



Number of thin bodies 

3 

NIB 

I 



Number of interference bodies 

4 

NSB 

I 



Number of slender bodies 

5 

NSUHF 

I 



Number of interpolation surfaces 

6 

NTBON 

I 



Number of thin bodies on the 
plane of sym. 

7 

NSBON 

I 



Number of slender bodies on the 
plane of sym. 

8 

NTBFON 

I 



Number of thin body elements on 
the P.O.S. 

9 

NSBON 

I 



Number of slender body elements 
on the P.O.S. 

10 

NTOTF 

I 



Number of force points total 
(NTBFP + NIBFP + NSBFP) 

11 

NTBFP 

I 



Number of thin body force points 

12 

NIBFP 

I 



Number of interference body 
force points 

13 

NSBFP 

I 



Number of slender body force 
points 

14 

NAFP 

I 



Number of aero force points 
(NTBFP + NSBFP) 

15 

NTOTC 

I 



Number of control points total 
(NTBCP + NIBCP + NSBCP) 

16 

NTBCP 

' I 



Number of thin body control points 

17 

NIBCP 

’ I 



Number of interference body 
control points 

18 

NSBCP 

I 



Number of slender body control 
points 

19 

NACP 

I 



Number of aero control points 
(NTBCP + NSBCP) 







Table 30. — (Continued) 









LABELED COMMON NAME; 


DESCRIPTION; FIT and buffer arrays for a file created 
with FETAD. 


ENG. 

NO. VARIABLE T DIMEN SION NOM. DESCRIPTION 

R 35 

R 514 








Table 30. — (Continued) 


LABELED COMMON NAME: CHKPRT 

DESCRIPTION! Checkout print option 

NO. 

VARIABLE 

H 

DIMENSION 

ENG. 

NOM. 

DESCRIPTION 

1 

1 

ICKPRT 

I 

1 

1 

1 

1 

Checkout print option 
=0 No checkout printing 
requested 

>0 Intermediate cheqkout print- 
ing was requested. 

ICKPRT is redefined for each 
primary overlay. 

In STRUCT(L217,2,0) ICKPRT 
triggers the following: 

*1 Print message upon entering 
and leaving each secondary 
overlay and print[M, ], 

[M2]/ and [M^] as generated. 

*3 Print items above plus in 
FORMAC(L221,2,2) print 
[CSDATA], {MCS}, [♦], [A], 

and [C]. 

Print items above plus a 
message from each VARDIM 
subroutine called. 
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Table 30. — (Continued) 
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LABELED COMMON NAME; 


CRDIMG 


DESCRIPTION: Storage array for last card image. 


NO, 


i VARIABLE 


T 


DIMENSION 


ENG. 

NOM. 


DESCRIPTION 


8 


Array to contain the last 80- 
column card image read. 


H 


Keyword extracted from ICARD. 
KEY contains the first 4 
characters of ICARD plus zero 
fill. 




Table 30. — (Continued) 
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Table 30. — (Continued) 


LABELED COMMON NAME: 


DESCRIPTION: Number of degrees of freedom 


NO. 


VARIABLE 


T 


DIMENSION 


ENG. 

NOM. 


DESCRIPTION 


1 


NDOF 


I 


Total number of degrees of 
freedom. 


2 


NDOFST 


I 


Number of structural degrees of 
freedom. 


3 


NDOFCS 


I 


Number of control surface 
degrees of freedom. 


4 NTDEF 


I. 


Number of mode used in DLPRES 
to develop gust forces at 
non-zero frequencies. 


5 


NPDEF 


I 


Number of mode used in DLPRES 
to develop gust forces at a 
frequency of zero. 
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Table 30. — (Continued) 


DUBFLS 

LABELED COMMON NAME: 

DESCRIPTION; Define files conta ining Doublet Lattice Data 


NO. 

VARIABLE 

T. 

DIMENSION 

ENG. 1 

NOM. 

DESCRIPTION 

1 

NGETP 

H 



Name of the Doublet Lattice 
geometry data file. 

2 

IGCASE 

I 



Geometry-^ case to be qsed. 

3 

NAETP 

H 



Name of the Doublet Lattice aero- 
dynamic data file. 

4 

I AC AS E 

I 



Number of the Doublet Lattice 
data case to be processed. 

5 

lAERO 

I 



Option from Doublet Lattice 






-1 only the [D] and [F] matrices 
are on NAETP 






=2 only the [Q] and [Cp] 
matrices are on NAETP 






=3 (D], [F], (Q], and (CP) are 

on NAETP. 
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EOMTPS 

LABELED COMMON NAME: 


DESCRIPTION: L217 (EOM) output file 


NO. 


1 

2 


VARIABLE 

EOMTAP 

EOMLOD 


!L 

1 

I 


DIMENSION 


W 2 










Table 30. — (Continued) 


LABELED COMMON 

FLXFLS 

NAME : 


DESCRIPTION: 

Define files containing 

FLEXSTAB data. 













NO. 

VARIABLE 

T 

DIMENSION 

ENG. 

NOM. 

DESCRIPTION 


GDTAPE 

IGCASE 

SDINDX 

SDDATA 

I AC AS E 

H 

I 

H 

H 

I 

1 


Name of the file containing 
FLEXSTAB geometry data. 

Number of data case on GDTAPE 
to be processed. 

Name of the file containing the 
index of the FLEXSTAB SDSSTP 
file. 

Name of the FLEXSTAB SDSSTP 
file, the file containing aero- 
dynamic data. 

Number of the data case on 
SDINDX to be processed. 
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LABELED COMMON NAME: 


GRDPEN 


DESCRIPTION: 


Define gradual penetration zones 









Table 30. — (Continued) 



The header array to be written 

onto 

the beginning of E0MTAP. 

Word 


Content 

1 

7HDYLOPLX 

2 

Program name & version 

3 

Date 

of run 

4 

NDOF 


5 

NLD. 


6 

NPAN 


7 

NFREQM 

8 

Q 


9 

V 


10-20 

Future use 

21 

Null 

(M^l 

22 

Null 

[M 2 ] 

23 

Null 

[M3] 

24 

Null 

[M 4 ] 

25 

Null 

[M5] 

26 

Null 

[Mg] 

27 

Null 

{C 2 > 

28 

Null 

{C,} 

29 

Null 

tftl 

30 

Null 


* The 

two arrays are equivalenced 










j 


Tab/e 30. — (Continued) 
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Table 30. — (Continued) 


LABELED COMMON NAME; 


MASSFL 


DESCRIPTION: Storage location of [M] 


NO. VARIABLE 


ENG. 

DIMENSION NOM. 


DESCRIPTION 


MTAPE 


MFILE 


MPRNT 


Name of. file containing [M], the 
generalized mass matrix. 

Number of logical fil'es to skip 
before reading [M]. 

Number of matrices to skip before 
reading [M]. 

Option requesting the printing 
of [M] if greater than 0. 
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Table 30. - (Continued) 


DESCRIPTION: Maximum Field Length 


DESCRIPTION 


Maximum field length used by 
the secondary overlay of a 
particular primary overlay. 

















Table 30. - 


LABE 

DESC 

LED COMMON 
RIPTION: 

NAME: OPTION 

List of program options 





NO. 

VARIABLE 

T_, 

DIMENSION 

ENG. 
NOM. _ 

1 

ISYM 

I 



2 

IDAMP 

I 



3 

HALFUL 

R 



4 

NGCP 




5 

NFREQM 

I 



6 

IQUASI 

I 



7 

I CENT 

I 




DESCRIPTION 


Key indicating symmetric 
(ISYM=0) or anti-symmetric 
(ISYM=1) analysis. 

Structural damping option 
(viscous damping, [M^], requested 
if IDAMP O) . 

HALFUL=1. indicates [M] and [K] 
are for full airplane. 

HALFUL=2- indicates [M] and [K] 
must be multiplied by 2. 

Number of gradual penetration 
control points (=1 means no 
gradual penetration) 

NOTE: (li=NGCPi=3 5) 

Number of frequencies at which 
aerodynamic data is defined 
(=0 means steady-state analysis 
is with KVALUE=0.) 

Key requesting Doublet Lattice 
quasi-steady AIC solution if >0 

Key indicating that FLEXSTAB 
control points are at thin body 
box centers if ICENT=3. At 
control points if ICENT=0. 




Table 30. — (Continued) 






Table 30. — (Continued) 
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Table 30. ~ (Continued) 








Table 30. — (Contmued) 


LAS PLED NAME: 


SCRFIL 







LABELED COMMON NAME: STIFFL 

DESCRI PTION : Storage location of [K] 



VARIABLE 


DntinuedJ 


_ DESCRIPTION 

Name of file containing [K], 

Number of logical files to skip 
before reading [K]. 

Number of matrices to skip 
before reading IF]. 

Option requesting the printing 
of (K} if greater than 0. 
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Table 30. - (Continued) 


LABELEf) COMMOM NAME: 


VARDIM 


DESCRIPTION: Accounting data for the VARDIM variable dimension 


array storage scheme. 


NO- 


VARIABLE 


T 


DIMENSION 


ENG, 

NOM. 


DESCRIPTION 


1 


NMAX 


I 


Maximum number of arrays to be 
defined with VARDIM. 


2 


NENTRY 


I 


Number of arrays currently 
defined- 


3 


LWAVAL 


I 


Last word available in blank 
common . 


4 


LWUSED 


I 


Last word used by VARDIM in blank 
common . 


5 


LKAT 


I 


6 


MAXUSD 


I 


Pointer to the VARDIM KATLOG 
array in blank common {=!). 

Maximum core address used to 
store a VARDIM array during the 
current program's (overlay's) 
execution . 
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Table 30. 


LABELED COMMON NAME: 


Blank 


DESCRIPTION: Variably dimensioned arrs 


NO. 


VARIABLE 


1 


1 BLANK 1 

* 


jlBLNK 




T 


DIMENSION 


ENG. 

NOM. 



J 


Storage 


DESCRIPTION 


The blank common array is loaded 
at the end of a program. The 
program uses the area between 
BLANK (1) and the end of field 
length for storage of variably 
dimensioned arrays - arrays with 
dimensions dependent upon the 
problem size. 


1 


4.0 EXTENT OF CHECKOUT 


Eadi module of L217 (EOM) was tested during development using numerical data cases 
which could be easily checked. The final version of the program was tested with a set of 
physically meaningful verification test cases. The verification cases and the options 
they exercise are displayed in table 31. 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
May 1977 


173 


Table 31. — Options Tested by Verification Data Cases 


Case 

Value 

or 

Option 

1 

FLEXSTAB 

N5KA 

2 

FLEXSTAB 

N5KA 

Anti-Sym 

3 

DUBFLX 

N5KA 


5 

DUBFLX 

N5KA 

6 

DUBFLX 

H7WC 

7 

DUBLFX 

N5KA 

Quasi- 

Steady 

8 

FLEXSTAB 

R7WC 

Document 

Sample 

Problem 

NDOFST 

5 

5 

j_ ... . 

5 

2 

5 

2 

5 

■ 

2 

NDOFCS 

0 

0 

0 

1 

0 

1 

0 

1 ■ 

SYM/ANTI 

SYM 

ANTI 

SYM 

SYM 

SYM 

SYM 

SYM 

SYM 

HALF/FULL 

HALF 

HALF 

HALF 

HALF 

HALF 

HALF 

HALF 

HALF 

[M] 

Cards 

Cards 

Cards 

Tape 

Cards 

Tape 

Cards 

Tape 

[K] 

Cards 

Cards 

Cares 

Tape 

Cards 

Tape 

Cards 

Tape 


Cards 

Cards 

Cards 

Cards 

Cards 

Cards 

Cards 

Cards 

C.S. Freedoms 

No 

No 

No 

T„,9 ,0 

No 


No 






z y . X 


Z' y .X 


Z y X 


FLXAIC 

FLXAIC 

DLPRES 

DLPRES 

DLAIC 

DLAIC 

DLAIC 

FLXAIC 

NTDEF/NPDEF 



1,2 

‘ 1,2 





IQUASI 





No 

No 

Yes 


ICENT 

No 

No 






No 

Gust Zones 

3 

3 

1 

1 

1 

1 . 

3 

1 

Print 

M,F,E* 

m,f,e 

M,F,E 

m,f,e 

M,F,E 

M,F,E 

M,F,E 

M,F,E 

Select Freq. 



1,2 

1,2 

1.2 

1,2. 



Slender Bodies 

2 

2 

3 

1 

3 

1 ■ 

3 

• 1 ' 

Elements 

12 

12 

9 

10 

9 

10 ; 

9 

10 

Thin Bodies 

3 

' 3 

3 

4 

3 

4 , 

3 

2 

Panels 

10 

,10 

10 

16 

10 

16 , 

10 

16 

Control Surface 

No 

No 

No 

No 

No 

No 

. No 

1 

Total Elements 

40 

40 

43 

62 

43 

62 

43 

. 62 

Modes 

5 

5 

5 

3 

5 

3 

5 

3 

Frequencies 

1 

1 

2 

2 

! 

2 

2 . 

2 

1 


* M,F,E means Modes, Forces, and Equations of Motion. 
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